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ABSTRACT 

Using the third data release of the Sloan Digital Sky Survey (SDSS) we have rigorously de- 
■ fined a volume limited sample of early-type galaxies in the redshift range 0.005 < z ^ 0.1. 

' We have defined the density of the local environment for each galaxy using a method which 

'nI" , takes account of the redshift bias introduced by survey boundaries if traditional methods are 

used. At luminosities greater than our absolute r-band magnitude cutoff of -20.45 the mean 
density of environment shows no trend with redshift. We calculate the Lick indices for the 
entire sample and correct for aperture effects and velocity dispersion in a model independent 
way. Although we find no dependence of redshift or luminosity with environment we do find 
[ that the mean velocity dispersion, a, of early-type galaxies in dense environments tends to 

be higher than in low density environments. Taking account of this effect we find that several 
indices show small but very significant trends with environment that are not the result of the 
correlation between indices and velocity dispersion. The statistical significance of the data is 
Q ■ sufficiently high to reveal that models accounting only for a-enhancement struggle to pro- 

|H I duce a consistent picture of age and metallicity of the sample galaxies, whereas a model that 

c/3 . also includes carbon enhancement fares much better We find that early-type galaxies in the 

field are younger than those in environments typical of clusters but that neither metallicity, 
a-enhancement nor carbon enhancement are influenced by the environment. The youngest 
early-type galaxies in both field and cluster environments are those with the lowest a. How- 
' ever, there is some evidence that the objects with the largest a are slightly younger, especially 

hj . in denser environments. 

Independent of environment both the metallicity and a-enhancement grow monotonically 
with (7. This suggests that the typical length of a the star formation episodes which formed the 
stars of early-type galaxies decreases with a. More massive galaxies were formed in faster 
bursts. 

We argue that the timing of the process of formation of early-type galaxies is determined 
by the environment while the details of the process of star formation which has built up the 
stellar mass are entirely regulated by the halo mass. These results suggest that the star for- 
mation took place after the mass assembly and favours an anti-hierarchical model. In such a 
model the majority of the mergers must take place before the bulk of the stars form. This can 
only happen if there exists an efficient feedback mechanism which inhibits the star formation 
in low mass halos and is progressively reduced as mergers increase the mass. 

Key words: galaxies: evolution - galaxies: elliptical and lenticular, cD - galaxies: abundances 



1 INTRODUCTION 

Present-day spheroids are thought to contain a large fraction, ~ 
60%-80%, of the total stellar mass in the Universe (Fukugita & 
Peebles 2004). They are thus fundamental tracers of the assembly 
of baryonic matter in the different environments generated from 
primordial density fluctuations. 



There are well established scaling relations indicating that the 
gravitational potential well is one of the key drivers in the for- 
mation and evolution of these galaxies. Suffice here to recall the 
colour-magnitude relation of early-type galaxies in clusters (Bower, 
Lucey & Ellis 1992), the relation between colour and central ve- 
locity dispersion and the fundamental plane (Bender, Burstein & 
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Faber 1993). In contrast, the role played by the environment is not 
yet clear 

The environment may affect galaxy properties in two distinct 
ways: first, the initial conditions from which galaxies form are dif- 
ferent and second their subsequent evolution may be affected by 
processes whose efficiency depends on environment. 

As for the initial conditions, the hierarchical clustering 
paradigm predicts a very simple dependence on environment, in 
the sense that objects in denser regions form earlier (Lacey & Cole 
1994). Subsequently the environment (whether we mean the local 
galaxy density or membership of a cluster) may affect galaxy prop- 
erties in several ways. 

In galaxy clusters, the large scale gravitational potential is 
directly responsible for so-called 'galaxy harassment' which de- 
scribes the cumulative effects of high speed encounters on the 
dynamics and star formation rate of the galaxies (Moore et al. 
1996). Furthermore, the intra-cluster medium can remove high 
filling factor ISM from disk galaxies (Gunn Sl Gott 1972). Ini- 
tial conditions and subsequent environmental effects are the most 
likely explanation for the existence of the colour magnitude rela- 
tion (Bower, Lucey and Ellis 1992), and the morphology density 
relation (Dressier, 1980, Mcintosh, Rix & Caldwell 2004; Balogh 
et al. 2004; Hogg et al. 2004; Croton et al. 2005), respectively. 

In the field however, the separation of the different aspects of 
the environmental influence is rendered much less clear by the ob- 
servational difficulty of obtaining homogeneous samples. Various 
authors, using samples of ~ 100 galaxies, find evidence that early- 
type galaxies in low density environments are younger and more 
metal rich than those in clusters (Longhetti et al. 2000, Thomas et 
al. 2005, Denicolo et al. 2005). Severeal authors (E.g. Kuntschner 
et al 2002, Denicolo et al. 2005, Thomas et al. 2005) also find that 
more massive and older galaxies have increased [a/Fe] values, sug- 
gesting more rapid star formation time-scales. Sanchez-Blazquez et 
al. (2003) use a sample of 98 galaxies to find lower C4668 and CN2 
index values in the centre of the Coma cluster than in the field. 

However, large and homogeneous galaxy samples are needed 
to provide the statistical discrimination necessary to identify small 
changes among stellar populations. Surveys such as the 2dF Galaxy 
Redshift Survey, Las Campanas Redshift Survey and the ongoing 
Sloan Digital Sky Survey (SDSS) have all been used for such stud- 
ies. A sample of 1823 galaxies (of all types) has been used by 
Balogh et al. (1999) to find stronger 4000 Angstrom break, and 
weaker [Oil] and H5 indices in the cluster environment. These au- 
thors conclude that the last star formation episode occurred more 
recently in low density environments. Bemardi et al. (1998), with a 
sample 931 early-type galaxies found only a minimal difference in 
the Mg2— (T relation, implying that field galaxies are only approxi- 
mately 1 Gyr younger than those in clusters. 

Bernardi et al. (2005a) have used a sample of over 39,000 
galaxies from the SDSS to find that age correlates with velocity 
dispersion from an analysis of the colour-magnitude-velocity dis- 
persion relation. Based on a similar sample Bernardi et al. (2006) 
analyse Lick indices to conclude that early-type galaxies in dense 
environments are less than 1 Gyr older than field objects and have 
a similar metallicity. Similar results have also been obtained by 
Smith et al. (2006) and Nelan et al. (2005). 

In this paper we study the absorption line indices of a large 
sample of early-type galaxies selected from the third release of the 
SDSS. We couple the high statistical significance of the sample 
with the most recent stellar population models to obtain an unbi- 
ased picture of the history of star formation in spheroids, as a func- 
tion of their velocity dispersion, in different environments. 



The paper is organised as follows. Section 2 is devoted to the 
sample selection and definition of environment. In Section 3 we 
derive the narrow band indices and describe the various corrections 
to transform them to the Lick IDS system. In Section 4 we analyse 
the data by means of our new simple stellar population models. 
Our findings are discussed in Section 5. Finally, in the conclusion 
Section we present an evolutionary scenario as a function of galaxy 
mass and environment. 



2 THE SAMPLE AND DATA 
2.1 Sample selection 

The primary sample used in this study was selected along similar 
lines to the sample used by Gomez et al. (2003) and Nikolic, CuUen 
& Alexander (2004), i.e., it is a luminosity- and (pseudo) volume- 
limited sample. The parent sample consisted of all objects, as of 
Data Release 3 (DR3), observed spectroscopically by the SDSS 
and classified as galaxies'. In practice, all of the galaxies that are 
part of our sample, given its selection criteria, should have been se- 
lected for spectroscopic observations as a part of the of SDSS Main 
Galaxy Sample (MGS) described by Strauss et al. (2002). Before 
selecting the volume-limited sample a small number of spurious en- 
tries were identified and removed: galaxies within 5 arc seconds of 
other galaxies (most likely de-blending errors), galaxies with red- 
shift confidence of less than 0.7 (they cannot be reliably placed in 
the volume limited sample) and galaxies with z-band magnitudes 
fainter than 22.8, i.e., not reliably detected in the z-band (these are 
most likely spurious detections). 

The volume limited sample was then constructed by selecting 
galaxies in the redshift range 0.005 < z ^ 0. 1 and with Petrosian 
absolute magnitudes brighter than -20.45. Since the r-band magni- 
tude limit of the MGS is 17.7, the absolute magnitude limit adopted 
should avoid a Malmquist bias in the volume limited sample. 

In order to limit our sample to include only early-type galaxies 
we then adopted the following selection procedure: 1) The com- 
pactness ratio as defined in the SDSS was constrained to be less 
than 0.33 (see Shimasaku et al. 2001). 2) Any object with detected 
emission in the Ha line was excluded in order to minimise the ef- 
fects of Balmer emission line infilling. This is particularly impor- 
tant for the correct interpretation of the Hj8 absorption feature. To 
this sample we further restricted the objects of study to, 3) those ob- 
jects with a well measured velocity dispersion, and 4) those objects 
where the density of the environment could be well determined, in- 
dependently of survey boundaries (see below). After such selection 
our sample consisted of 3614 galaxies. In section 4.; we describe 



a further selection based on the signal-to-noise ratio of the specttra 
to be included in the final regression analysis. 



2.2 Estimating density of environment 

The algorithm used to quantify the density of the environment is 
similar to that used by Dressier (1980) in his original study of the 
relationship between density and morphology: it is based on find- 
ing the projected separation to the n-th nearest neighbour, r„. This 
distance can, if required, be used to estimate the surface density of 
galaxies via Ob „ =njr^n. 

Unlike the original algorithm, the one presently employed 



' These ar e all of the objects in files gals-DR2.fit.gz and gals-DR3.fit.gz 
available at http://das.sdss.org/DR3/data/spectro/ss_tar_23l 
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makes use of the measured redshifts of galaxies to search for near- 
est neighbours only in the redshift shell Av = ± 1200km s~' around 
the galaxy being considered following Gomez et al. (2003); Goto 
et al. (2003) (although our maximum velocity difference is slightly 
larger). Because of this modification, no correction for the back- 
ground galaxy field density need be made. Since in our approach 
redshift information is required, the search for «-th nearest neigh- 
bour is necessarily restricted to a sample of galaxies with measured 
redshifts. To avoid a dependence of the calculated density scale on 
redshift, we search for n-th nearest neighbour in the volume-limited 
catalogue itself. A photometric approach to calculating galaxy den- 
sity in the SDSS is discussed by Hogg et al. (2003). 

Since density of environment is a non-local measure, survey 
boundaries will inevitably present a problem. Galaxies close to 
one of the boundaries will have the separation to their «-th neigh- 
bour (the density metric that we are using in this study) overesti- 
mated and therefore their density of environment underestimated; 
this effect is quite clearly illustrated in the simulations presented by 
Miller et al. (2003). Hence galaxies close to survey edges have un- 
reliably measured densities and should not be considered in further 
analysis. The non-compact nature of the SDSS survey area makes 
this effect quite pronounced, although less so in DR3 than in earlier 
releases. 

The usual approach to dealing with this problem is to remove 
from further analysis galaxies with n-th neighbours at a larger sep- 
aration than a survey boundary. This was shown by Miller et al. 
(2003) to be effective, but such an approach also preferentially re- 
moves galaxies in low density environments at low redshift com- 
pared to high redshift, as shown in the upper panel of Fig. 0. This 
can be understood by considering a population of galaxies with a 
fixed physical distance to their n-th neighbour: with increasing red- 
shift, the angular separation to a survey boundary can be smaller 
and still allow good estimation of density; therefore the necessary 
margin around the survey boundary decreases and the effective use- 
ful area of the survey increases with increasing redshift. 

The correlation between density and redshift in the sample af- 
ter such an edge-correction procedure can be avoided by removing 
a somewhat larger number of galaxies as follows. If r„ is the sep- 
aration to the n-th neighbour and is the separation to a survey 
edge, only galaxies which satisfy r„f^ ^ should be retained in 
the edge-corrected sample, where zq is a fiducial redshift which 
should also be (close to) the lower redshift limit of the sample. 

Finally, the choice of n to use as a density metric will also de- 
termine how much of an effect survey boundaries have. We follow 
Goto et al.(2003) and use n = 5 throughout. The edge-corrected 
sample is shown in the lower panel of Fig. |lj 

Fig. ^ shows the distribution of our measure of environment 
for all the galaxies in our sample. 

In Fig. ^ we compare our measure of the local galaxy density 
with the number of cluster members for those Abell clusters where 
our sample contains at least one galaxy. There is a correlation be- 
tween the measures but one is not a good predictor of the other, and 
neither would one especially expect them to be. Our measure of 
environment quantifies the local environment of a galaxy whereas 
the total number of galaxies in an Abell cluster is a measure of 
the large-scale environment (recall also that we define environment 
using only those galaxies bright enough to be detected at all red- 
shifts in our sample). The most important point about Fig. ^ is that 
it shows that galaxies in our sample that have a density parameter 
above 1 are in environments typical of galaxy clusters. 38% of our 
galaxies have a value of 1 /r^ > 1 .0, and 10% l/r^ > 2.0. 

As a further check for the reliability of our measure of envi- 




Figure 1. The distributions — represented as the median (solid line) and the 
25 and 75 percent quartiles (dashed lines) — of density as a function of red- 
shift for samples corrected for the effects of survey boundaries in the tradi- 
tional way (top) and using the technique described in the text (bottom). 
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Figure 2. The distribution of the environmental density for our sample. 



ronment we visually checked the fields of our sample galaxies with 
the 20 highest and 20 lowest values of l/r^- This confirmed that 
low values do indeed represent galaxies in poor environments and 
visa-versa. 



3 CALCULATING THE LINE INDICES 

Since one of the final goals of this study is to interpret the indices 
with models based on the Lick system we first transform the SDSS 
spectra to the wavelength dependent resolution of the Lick-IDS 
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Figure 3. Comparison of our measure of environment, l/fs, witli tlie num- 
ber of cluster members as given in Abell. The values of l/rs ^re the mean 
value for the galaxies in our sample associated with each cluster. The Abell 
cluster designation is shown near each point. 



spectra. This corresponds to ~ 11.5 A at 4000 A, ~ 8.4 A at 5000 A 
and ~ 9.8 A at 6000 A. 

From the SDSS spectra we then extract the 21 line-strength 
indices of the original Lick-IDS system (Trager et al. 1998) plus 
the additional indices HyA, H5A and HyF, H5F (Worthey & Otta- 
viani 1997) and the indices B4000 (Hamilton, 1985) and HK (Rose, 
1984, 1985). Our extraction makes use of the redshift of each ob- 
ject as provided in the SDSS. For this extended Lick system we 
transform the passband definition of Trager et al. (1998) to the vac- 
uum (SDSS spectra are given in terms of vacuum wavelengths). A 
comparison with the indices in common with the associated SDSS 
tables show that they lie on a one-to-one relation. 

In our analysis the transformation to the Lick system, which is 
usually made by comparing selected observed standard stars with 
the Lick IDS observations, is not possible. This is simply because 
of the lack of observations of standard Lick stars in the adopted 
SDSS data release. For this reason we work with index differences, 
i.e. we do not consider absolute values but only trends with respect 
to a reference population. This strategy avoids significant biases in 
our analysis, because it has been shown that the transformations to 
the Lick IDS standard stars mostly consist of an offset correction 
(see e.g. Puzia et al. 2002, Rampazzo et al. 2005). Our results will 
also be much less dependent on the zero point calibration of the 
SSP models. 

The errors on each index measurement were calculated by first 
measuring the rms noise in each spectrum and then reproducing a 
model spectrum with the same rms. We then used a Monte-Carlo 
method to measure the indices for 200 realisations of the spectrum. 
We define the error on our measurement of an index as the standard 
deviation of the index values so obtained. 



3.1 Velocity dispersion correction 

Because we are interested in the intrinsic index- CT relation for early- 
type galaxies it is important that no bias be introduced by an inap- 
propriate velocity dispersion correction. However, as we show in 
the appendix, although most indices show some correlation with 
the velocity dispersion of a galaxy, and velocity dispersion is in 
turn correlated with environment, our results on the influence of 
the environment are not sensitive to our velocity dispersion correc- 
tion. The velocity dispersion correction is thus important only for 



the shape of the index-C7 relations. The velocity dispersion for each 
galaxy is provided in the SDSS. 

To ensure that our derived index-a relations are not a result 
of an error in velocity dispersion correction we have made the cor- 
rection in two different ways and carried out our analysis with both 
versions. 

In the first we simply smooth all the spectra to a velocity res- 
olution of 400 kms^'. Although this will result in index values 
which are all lower than they would be at say the resolution of the 
Lick system, we are sure that we have not introduced any model 
dependent effects. 

In the second method we have smoothed the spectra of G and 
K-type giant stars observed by the SDSS to a number of different 
velocity dispersions from zero to 400 kms^' . For each index of a 
given galaxy we identify a stellar spectrum with the same velocity 
dispersion and with similar line index values, and then use this star 
to make the correction to the Lick IDS system. 

A comparison between these 2 methods has been made by us- 
ing the stars observed by the SDSS to correct indices to a velocity 
dispersion of 400 kms^' and comparing to the results of the first 
method. Only in some of the bluer indices (Ca4227, G4300, Fe4383 
and Ca4455) was any difference seen in the resulting index-CT rela- 
tion. 



3.2 Aperture correction 

Every spectra in the SDSS is taken through a fibre with a diameter 
of 3" and so the physical radius sampled for any given galaxy is a 
function of distance. We must apply an aperture correction to re- 
move this effect. Rampazzo et al. (2005) have measured the spatial 
index gradients in 50 elliptical galaxies and we have used the mean 
radial gradients found for each index to correct our index values to 
the standard radius of re/10 where re is the effective radius. Using 
the data of Rampazzo et al. we take {I{re/2) -I {re/ 10)) /I (re/ 10) 
as a measure of the spatial index gradient and plot this against log CT 
for all 50 galaxies in their sample. A least-squares, straight line fit 
to these plots is then found; 



(/(r,/2)-/(r,/10)) 



:xlog(a/100)+>' 



(1) 



Hre/lO) 

The values of x and y for each index are reported in table |l| The 
SDSS gives values for the half light radius (re) of each galaxy and 
so we can correct the index measured with a 3" SDSS fibre to the 
standard radius of (rg/lO) by assuming a linear change in index 
value with radius; 



I{re/W) 



/(SDSS) 



(xlog(cT/100)+y) 



+ 1 



(2) 



where /(SDSS), the index measured from the SDSS spectrum, is 
the index measured within a radius of l"5. CT is in units of kms^^ 
and rg in arcseconds. 

Fig. ^ shows the aperture corrections applied for each index. 
The accuracy of the aperture corrections depends on the accuracy 
of the value for the half light radius reported in the SDSS. This in 
turn will depend principally on the seeing. The fractional error will 
therefore increase with redshift, on average, as the angular sizes 
of the galaxies decreases, with the radii tending to be increasingly 
overestimated at greater distances. This effect has yet to be quan- 
tified in the current release of the SDSS and the effect would in- 
fluence the gradients of our index- a plots (Fig. However, as we 
see below, our derived index- C7 relations agree very well with those 
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Table 1. Parameter values for fits to the radial index variations measured by 
Rampazzo et al. (2005). Values for x and y in equation |l| are given. 



Index 


X 


y 


Index 


X 


y 


CNl 


0.468 


-0.434 


Fe5270 


0.195 


-0.0937 


CN2 


0.658 


-0.463 


Fe5335 


0.0340 


-0.0565 


Ca4227 


0.441 


-0.0999 


Fe5406 


-0.0990 


-0.0204 


G4300 


0.129 


-0.0531 


Fe5709 


0.131 


-0.0225 


Fe4383 


-0.0344 


-0.0570 


Fe5782 


0.227 


-0.0986 


Ca4455 


-0.111 


-0.0716 


NaD 


0.127 


-0.193 


Fe4531 


0.103 


-0.0478 


TiOl 


0.143 


-0.108 


C4668 


0.318 


-0.219 


Ti02 


0.0434 


-0.0573 


Hi3 


0.410 


-0.0369 


B4000 


-0.0785 


0.0574 


Fe5015 


0.517 


-0.158 


HK 


0.0165 


0.00396 


Mgl 


0.132 


-0.142 


HDF 


1.18 


-0.539 


Mg2 


0.131 


-0.107 


HGF 


0.126 


-0.183 


Mgb 


-0.0068 


-0.0594 









derived from other studies based on more local samples and so we 
expect the influence of a seeing-induced effect to be small. Further, 
we note that because our various sub-samples have similar redshift 
distributions (Fig. ^ our results on the influence of environment are 
insensitive to the gradient of the index-(j relations (see appendix). 



4 RESULTS 

Below we quantify our results making use of both simulations and 
models. We stress however that the results are qualitatively appar- 
ent in the data before the use of either. In the logcT-index plane, 
for example, several indices show a similar gradient for different 
environments but with an offset in index value. This immediately 
implies an effect of environment on the index that is not due to the 
elevated a in denser environments (see section 4.2). Because we 



have taken great care to minimise selection bias we are able to treat 
the sample as a whole for the statistical analysis. 

The velocity dispersion, a, does not have a one-to-one rela- 
tion to the mass of the galaxy. Recent work by Cappellari et al. 
(2006) has investigated the relation between galaxy mass-to-light 
(M/L) ratio and the line-of-sight component of the velocity dis- 
persion within the effective radius. They provide relations (their 
equations 7 and 10) which allow a transformation from measured 
velocity dispersion to galaxy mass; 



Mio = (16.5 ±7.8) CT; 



3.11±0.43 
200 



(3) 



where Miq is the galaxy mass in units of 10'"M,r; and 02(10 is the- 
luminosity weighted mean velocity dispersion within one effective 
radius in units of 200 km s^ ' . The errors have been propagated from 
the parameter values found by Cappellari et al. (2006). We do not 
make any conversion to mass but those who wish to can apply this 
formula. 



4.1 Trends of line indices with velocity dispersion 

Figure ^ shows the index-CT relations for our entire sample after 
the aperture and velocity dispersion corrections have been made. 
The gradients of the relations are reported for a simple first order 
fit along with the corresponding errors. Most indices show some 
correlation with the velocity dispersion. 

The following indices show a significant (> 5c7) posi- 
tive correlation: CNl, CN2, G4300, Fe4383, Ca4455, Fe4531, 



C4668, Mgl, Mg2, Mgb, Fe5335, Fe5406, NaD, TiOl, Ti02, 
<Fe>, [MgFe] (Gonzalez, 1993) and [Mg/Fe] (defined as 
Mgb/{Q.5{Fe5210 + Fe5335))). Negative gradients are instead 
seen for: Hj8, H5, Hy, Fe5015, B4000 (note the definition of this 
index in the SDSS is inverted with respect to that of many authors) 
and HK. No appreciable gradient is evident for: Ca4227, Fe5270, 
Fe5709 and Fe5782. Below we consider straight line fits of the form 
/ = ao + 01 log a or log/ = Aq +A\ logCT. 

Most work on the trends of indices with velocity dispersion 
has been done for the Mg2-(7 relation for which we find a\ = 
0.209 ±0.0056 and ao = -0. 198 ± 0.013. Bernardi et al. (1998) 
find ai = 0.224 ± 0.008 and ao = -0.23 ± 0.019 for a sample of 
931 galaxies and Kuntschner et al. (2002) find corresponding val- 
ues of 0.216 ±0.025 and -0.236 ±0.057 for Fornax spheroidals. 
Both studies give fits consistent with ours. Jorgensen (1997) finds 
values of a\ = 0.196 ±0.009 and ay = —0.155. The slope is in 
agreement with our value whereas the intercept seems to be slightly 
lower (although this author quotes no error on the intercept value). 

For MgB we find fli = 3.78 ±0.153 (Ai = 0.389 ± 0.0165) 
which agrees well with Bernardi et al. (2003) and is similar to the 
relation found by Thomas et al. (2004) (although they do not ex- 
plicitly quote the gradient). 

For <Fe> we find ai =0.449 ± 0.072 and qq = 1.78 ±0.162 
(Al = 0.0820 ±0.0118 and Aq = 0.258 ± 0.0267) respectively. 
These values are both consistent with that found by Bernardi et 
al. (2003) and Jorgensen (1997). 

Our [Mg/Fe]-a relation has ai = 1.04 ± 0.0700, ao = 
-0.747 ±0.159 (Al = 0.306 ± 0.0194, Aq = -0.499 ± 0.044). 
Kuntschner et al. (2002) find values in the range A 1 --^0.18-0.30 
for < 100 early-type galaxies in Virgo, Coma and Fornax. Our 
value is consistent with their findings despite the fact that our 
value refers to early-type galaxies in a wide range of environments. 
Our [Mg/Fe] -CT relation is marginally consistent with Bernardi et 
al. (2003) who find Ai = 0.22 ±0.04. Trager et al. (2000) find 
Al = 0.33 ± 0.01 and Aq = -0.58 ± 0.01 which is similar to the 
gradient we find but with a slightly lower intercept. 

We obtain a slope for the Hj8 index of aj = -1.75 ±0.0670 
(Al = -0.450 ±0.0191). This is similar to that found by Thomas 
et al. (2005) (~ —1.7) but is somewhat steeper than that found by 
Bernardi et al. (2003), who obtain d(logH/3)/d(log cr) = -0.24 ± 
0.03. 

A thorough treatment of the index-CT relations in cluster galax- 
ies has been carried out by Nelan et al. (2005) (their table 7). Most 
of the gradients measured by these authors are similar to the gradi- 
ents we find for our high density sub-sample, but some indices show 
large differences from our measured trends. These are Ca4227, 
C4668, Fe5015 and Fe5270. We note that apart from the differences 
in sample characteristics these authors apply an aperture correction 
to a fixed physical radius that is not a function of ct; our correction 
instead is to re/ 10 and is applied as a function of velocity disper- 
sion. Aperture corrections have not traditionally taken into account 
the fact that radial index gradients are a function of CT. We note that 
the indices in which we see the most significant disagreement tend 
to have aperture corrections that vary more strongly with CT (Fig.^. 

In section 4.3, these relationships are interpreted in terms of 



physical parameters of the galaxies. 



4.2 The effect of environment on the line indices 

Having corrected our catalogue for the various biases discussed 
above we then divided it in terms of the density parameter 1 /rg to 
investigate the effect of environment on the line indices. Of course. 
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Figure 4. Aperture corrections as a function of velocity dispersion, a. The points are the data, and the solid and dashed lines show the correction applied for 
an effective radius of 2" and 4" respectively. 
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Figure 5. Index-C7 relations. The data have been corrected for both aperture and velocity dispension effects. The black line shows a least squares, straight line 
fit to the whole sample, the parameters of which are reported at the top of each panel. The red and green hues show fits to low and high density sub-samples 
with l/r'i < 0.5 and l/r'; > 1.5 resoectivelv. 
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Figure 7. Mean value of the velocity dispersion, (J, as a function of our 
measure of the environmental density, l/r^. The error bars are the standard 
deviation of the distribution of velocity dispersion divided by the square 
root of the number of objects. 



if the distribution of our sub-samples differs in any parameter that 
is correlated with the index we must take account of this. 

We divided our total sample of 3614 galaxies in the follow- 
ing way. First we selected a low density sample made up of those 
galaxies with \/r^ < 0.5 (726 objects) against which to compare 
galaxies in increasingly dense environments. These were then de- 
fined as follows: 1) 0.5 < l/rj < 1.0 (1513 objects), 2) l/rj > 0.5 
(2888 objects), 3) 1 /rj > 1.0 (1362 objects), 4) 1 /rj > 1.5 (706 ob- 
jects), and 5) l/rj > 2.0 (357 objects). The upper density bound on 
the first of these was intended to search for effects at environmental 
densities below those typical of clusters. 

In Fig. P we show the distributions of velocity dispersion, red- 
shift and luminosity for each of the sub-samples compared with the 
sub-sample with \/r^ < 0.5 along with the probability that the two 
distributions are drawn from the same parent population. We note 
that the luminosity and redshift distributions are similar in each 
case (although the two lowest density bins for the redshift distribu- 
tion have a probability of only 0.7% of being drawn from the same 
parent population their means are similar). The distributions of ve- 
locity dispersion on the other hand are very significantly different in 
all cases. The velocity dispersion tends to be higher in denser envi- 
ronments and increases monotonically with our measure of the den- 
sity of environment. At l/rs, < 0.5, <a>= 182.7 kms^', whereas 
for l/rs > 2.0, <C7>= 193.9 kms"'. This trend is illustrated in 
Fig. In order to test that the differences in velocity dispersion 
distributions are not due to a special distribution in the sub-sample 
with 1 /rg < 0.5 we also compared each environment with the sub- 
sample with 0.5 < 1/^5 < 1.0. Significant differences remained (al- 
though at lower significance due to the smaller total number of ob- 
jects and the smaller differences in 1 /rg. 

As we saw above, most indices show some correlation with the 
velocity dispersion. The differences in mean velocity dispersion for 
different environments therefore mean that differences in the index 
values between environments will have a contribution due only to 
the index-CT relation. The two effects might entirely explain any 
index differences we find with environment. In order to take into 
account this effect we have employed a Monte-Carlo approach. 

Using the entire sample (before defining regions in different 
environments) we consider the index-CT relation. Binning the data 
in CT with a bin width of 10 kms^' we plot the value of the mean 



index per bin against CT and linearly interpolate between points. 
This line is used to define the index- CT relation. 

For a given sub-sample of environmental density we then take 
the measured velocity dispersion value of each galaxy and translate 
it into a value for the index via the index-c relation. In the process 
we alter the derived index value by a random amount drawn from 
an appropriate Gaussian distribution so that the observed standard 
deviation of measured indices is correctly reproduced. We thus cre- 
ate a model data set with an index distribution defined only by the 
velocity dispersion of the sub-sample and the index-c relation. 

In this way we generate a model index distribution for a given 
sub-sample of environmental density. Any difference between the 
index values from two model sub-samples defined in this way is due 
only to the difference in the distribution of velocity dispersion. We 
quantify this difference by subtracting the means of two environ- 
mental density sub-samples for 1000 realisations of the model. The 
distribution of this model difference gives us a distribution against 
which to test the observed difference in index. If the real data show 
a difference well in excess of that seen in the simulated data then 
we identify an effect of environment on the index. 

Fig. |8| displays the normalised index differences for 5 environ- 
ments of increasing density, each compared to the sub-sample with 
the lowest density (l/fg < 0.5). The normalised index difference is 
given by. 



A < 7.5 



</.!>-< /o > 
</ > 



(4) 



where < /j > is the mean value of the index in the sub-sample with 
mean environmental density s, Iq is the mean value of the index in 
the comparison sub-sample and < / > is the mean index value for 
the whole sample over all environmental densities. Table ^ sum- 
marises the results of Fig ^. 

The following indices show a statistically significant density 
dependence in the sense that index values are higher in more dense 
environments: CNl, CN2, Ca4227, G4300, Fe4383, (Ca4455), 
Mgl, Mg2, Mgb, (Fe5335) NaD, Ti02, and [MgFe] (parentheses 
indicate that we find an offset but not a trend with density). The 
B4000 index shows the opposite trend with lower values in more 
rich environments but this is simply due to how the index is defined 
in the SDSS. Although none of the differences attributable to en- 
vironment significantly increase the scatter of index values of the 
population as a whole, the differences are very significant. 

Indices which show an offset but no trend with density over 
the range of the plots of Fig. ^ show an excess of low values in the 
comparison low density sub-sample. 

The above indices have different sensitivities to age and metal- 
licity variations and, in general, an increase of a metal index may 
indicate increased age and/or larger metallicity. As for the effects of 
the metallicity, we recall that the sensitivity to separate elemental 
variations is also important (enhancement) and has been investi- 
gated by various authors (Tripicco & Bell 1995, Trager et al, 1998, 
Thomas, Maraston & Bender, 2003, Korn, Maraston & Thomas 
2005, Annibali et al. 2005). Indeed, from the indices in which an 
effect is seen with density of environment (and from those which 
show little or no effect) we can also attempt to identify elemental 
abundance variations. In this respect carbon is of particular impor- 
tance. 

Tripicco & Bell (1995) and Korn et al. (2005) showed that 
carbon has a strong influence on many of the indices in which we 
find a trend. However it is striking that C4668, in spite of showing 
the greatest sensitivity to carbon, displays at most only a marginal 
increase with environmental density. This is inconsistent with an 
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Figure 6. Distribution of velocity dispersion, redshift and z-band luminosity for various density environments. In each case the dotted hne is the distribution 
for the comparison sub-sample with l/cs < 0.5 and the solid (red) line is that for increasingly dense environments. Density increases from top to bottom in 
each panel. The density of the environment considered, number of galaxies, and the KS- probabihty that the 2 distributions are drawn from the same parent 
population are shown in each plot. Left: Velocity dispersion, middle: redshift, right: luminosity. The clear increase of the velocity dispersion with density of 
environment is quantified in Fig. ^ 



increase in carbon abundance being the cause of the differences we 
see for these indices. It also suggests that, if a general abundance in- 
crease is the cause of the index difference, carbon is a bad tracer of 
overall abundance and that its relative abundance should decrease 
as metallicity increases. 

We note that we do not reproduce the results of Sanchez- 
Blazquez et al. (2003) who find both a lower C4668 and CN2 index 
in for early-type galaxies in the central regions of the Coma clus- 
ter when compared to the field. This study was based only on Coma 
galaxies and the mean environmental density of galaxies in our sub- 
sample with l/rs ^ 1 .5 is unlikely to be as high as for those at the 
centre of the Coma cluster. 

Besides the various dependencies of indices on elemental 
abundances they are also, of course, affected by age. In order to 
make sense of the trends of the indices we must make use of the 
sensitivity to the different physical parameters simultaneously. 

4.3 Age, metallicity and a -enhancement 

In order to extract trends in age, metallicity and a-enhancement 
from our measured trends in the line indices we have used the re- 
cent models computed by Annibali (PhD thesis and Annibali et al. 
2005) and a multiple-linear regression analysis. 

These models account for the effects of age, metallicity and 
a-enhancement variations. The models have been computed by 
adopting the fitting functions that describe the dependence of the 
Lick indices on stellar parameters. The effects of the element par- 
tition has been introduced following the sensitivity to different el- 



ements provided by Kom et al. (2005). In computing the effects of 
a-enhancement the elements N, O, Ne, Na, Mg, Si, S, Ca, Ti were 
assigned to the enhanced group while Cr, Mn, Fe, Co, Ni, Cu, Zn 
to the depressed group. Because our analysis makes use of indices 
that show a mild to strong dependence on the relative abundance 
of carbon, we have also considered the effects of enhancing this 
element alone. 

Our selection of indices to include in the regression analysis is 
restricted to those satisfying two criteria. Firstly, the indices must 
be well modelled by the SSPs of Annibali et al. (2005) and sec- 
ondly, those that have fractional errors < 20%. In common with 
similar models by other workers (e.g. Thomas, Maraston & Bender 
2003), certain indices are not reproduced well by the adopted SSP 
models. The CN indices are underpredicted compared with globu- 
lar cluster data unless nitrogen is increased by a factor 3 with re- 
spect to the a-elements. NaD is sensitive to interstellar absorption 
(Burstein et al. 1984) and so is of limited use for stellar population 
studies (Worthey et al. 1994). The model predictions of the Ca4227 
index are slightly too high while Ca4455 suffers from calibration 
problems (Maraston et al. 2003). The TiOl and Ti02 indices ap- 
pear to be poorly calibrated, and furthermore measured TiOI in- 
dices may be affected by the presence of the broad NaD feature 
(Denicolo et al. 2005). Among the Fe indices, we adopted the clas- 
sic Fe5270 and Fe5335 indices plus the blue Fe4383 and Fe4531 
indces which are both sensitive to iron and are well calibrated.The 
other iron indices are relatively poorly calibrated (Fe5782), are less 
sensitive to Fe (Fe5709), or may be affected by emission contam- 
ination (Fe50I5, Kuntschner et al. 2002). For the age-sensitive in- 
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Figure 8. Comparison of normalised index values with the comparison low-density sub-sample (l/rs < 0.5) for increasingly dense environments. The hor- 
izontal axis shows values of the mean density of the environment ( 1 /rs ) and the vertical axis shows the fractional deviation of the index value with respect 
to the sample mean. Vertical en'or bars show the results of Monte-Carlo simulations which take into account the index— CT relation and the difference in the 
distribution of velocity distribution between environments of different density. Open diamonds represent the observed index difference. For each index four 
panels are shown. The left panel shows results for the whole sample. The other 3 panels show the results for 3 ranges of velocity dispersion. From left to right: 
a < 160 kms"', 160 ^ fj ^ 230 kms"' and a > 230 kms"'. 
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Table 2. Comparison of normalised index values for the highest and lowest 
density sub-samples shown in the left panel of Fig. ^ The data and the 
simulation refer to the whole sample. LDE and HDE refer to densities of 
1 /rj =0.712 and 1 /rj = 2.734. 



Index 



Simulation 
LDE 



Data 
LDE 



Simulation 
HDE 



Data 
HDE 



CNl 

CN2 

Ca4227 

G4300 

Fe4383 

Ca4455 

Fe4531 

C4668 

H/5 

Fe5015 

Mgl 

Mg2 

Mgb 

Fe5270 

Fe5335 

Fe5406 

Fe5709 

Fe5782 

NaD 

TiOl 

Ti02 

B4000 

HK 

iFe6 

[MgFe] 
[Mg/Fe] 
HdF 
HgF 



0.034±0.020 
0.015±0.013 
0.000±0.007 
0.000±0.003 
0.003±0.()05 
0.008±0.008 
0.003±0.006 
0.0()2±0.()05 
-0.009±0.008 
-0.003±0.005 
0.014±0.()07 
0.007±0.004 
0.010±0.007 
-0.001±0.005 
0.003±0.()06 
0.0()4±0.()08 
-0.006±0.012 
0.007±0.014 
0.011 ±0.006 
O.OlOiO.OlO 
0.007±0.005 
-0.001±0.002 
-0.001 ±0.002 
0.001±0.006 
0.006±0.004 
0.009±0.008 
-0.048±0.046 
-0.019±0.015 



0.066±0.028 
0.028±0.018 
0.020±0.010 
0.003±0.005 
0.017±0.007 
0.041±0.011 
0.018±0.008 
0.013±0.007 
-0.014±0.012 
-0.008±0.007 
0.032±0.010 
0.017±0.006 
0.033±0.009 
-0.001±0.008 
0.019±0.009 
-0.014±0.011 
-0.025±0.017 
0.003±0.020 
0.020±0.009 
0.027±0.014 
0.032±0.007 
-0.006±0.003 
-0.004±0.002 
0.009±0.008 
0.019±0.006 
0.027±0.011 
-0.237±0.066 
-0.062±0.022 



0.083±0.026 
0.039±0.017 
0.006±0.010 
0.000±0.004 
0.012±0.007 
0.015±0.012 
0.006±0.009 
0.009±0.007 
-0.025±0.011 
-0.007±0.007 
0.031±0.0100 
O.019±O.OO5 
0.023±0.008 
-0.001±0.008 
0.009±0.008 
0.010±0.012 
-0.011±0.015 
0.014±0.017 
0.026±0.009 
0.027±0.014 
0.018±0.007 
-0.004±0.002 
-0.001±0.002 
0.004±0.008 
0.014±0.006 
0.018±0.010 
-0.152±0.054 
-0.042±0.019 



0.227±0.036 
0.131±0.024 
0.047±0.015 
0.027±0.007 
0.044±0.010 
0.048±0.016 
0.008±0.012 
0.019±0.010 
-0.041±0.015 
0.002±0.010 
0.082±0.014 
0.048±0.007 
0.066±0.012 
0.016±0.012 
0.027±0.012 
-0.005 ±0.017 
-0.050±0.021 
-0.020±0.025 
0.049±0.012 
0.056±0.019 
0.064±0.009 
-0.017±0.003 
-0.006±0.003 
0.022±0.012 
0.044±0.008 
0.044±0.015 
-0.480±0.082 
-0.157±0.027 



dicators we adopted the classic H/3 index, which is mostly inde- 
pendent of a-enhancement effects but may be affected by emission 
contamination and we preferred H5 to Hy as it is less affected by 
nebular emission. 

Our final list of indices is: H/3, H5, Mgl, Mg2, Mgb, Fe4383, 
Fe4531, Fe5270, Fe5335, G4300 and C4668. (We note that the in- 
clusion of Hj8 has no significant effect on the results. The primary 
reason for it's inclusion was to allow comparison with the results of 
other authors.) We find that our results are also not a strong function 
of the adopted S/N cutoff of the spectra. 

We express the observed index variation as a linear combina- 
tion of age, metallicity, a-enhancement and carbon variation. 



5/: 



dl 



(9log(f 
dl 



■8\og{t) 



dl 



■5[a/Fe] 



^log(Z) 
dl 



5log(Z)- 
5[C/H] 



(5) 



a[c/H] 

where the partial derivatives are derived by the models and the ob- 
served variation dl is obtained with respect to a reference point. 
Our use of a linear combination of model parameters is justified 
in Fig. ^ where we see that the index variations are indeed linear 
with respect to changes in the model parameters in equation ^For 
a given observed data point, by considering any number of indices 
larger than four at once, we can perform a multi-regression anal- 
ysis and obtain the corresponding variations of 5Log(t), 5Log(Z), 
5[a/Fe] and 5[C/H]. For the purposes of the regression analysis 
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Figure 11. The number of objects in each bin of Fig. |l^ The symbols have 
the same meanings as in Fig. id. 



we have first considered the whole sample regardless of environ- 
ment and then 2 sub-samples chosen to be representative of the 
'field' and 'cluster' environment. For the field sub-sample we take 
1/^5 ^ 0.5 and for the cluster sub-sample we take 1/^5 ^ 1.5. We 
then bin the index data in a. Index variations are computed relative 
to the index values for the whole sample at C7 = 200 kms^' . 

In a preliminary analysis we have not considered the varia- 
tion of the abundance of carbon [C/H]. This allows a comparison 
of the results obtained with our new models with those obtained 
with the models of Thomas et al. (2004). The outcome of this anal- 
ysis can be summarised as follows. Although we always find an 
a-enhancement that increases with velocity dispersion our results 
for the age and metallicity depend on which indices are included in 
the analysis. The results are similar if we make use of the models of 
Thomas et al. (2004). Specifically, our results change with the in- 
clusion/exclusion of indices which have a strong dependence on the 
carbon abundance such as C4668 and Mgl. As we noted above the 
lack of strong C4668-(7 relation suggests that the index does not 
trace overall metallicity any more than the mean iron abundance. 
The different sensitivities specifically to carbon abundance of many 
of the indices makes the consideration of this element of particular 
importance. We therefore subsequently modified our analysis to in- 
clude independently the effects of carbon enhancement. Although 
the correction for carbon abundance is really just the most impor- 
tant part of a correction that should take account of other elements 
such as titanium and sodium we consider here the refinement only 
for carbon. As a consequence of the modification, our results for 
the age and metallicity using the refined model are much less de- 
pendent on the indices included in the regression analysis. With 
this modified model we can include a larger number of indices and 
exclude only those which have a relatively strong dependence on 
elements not explicitly modelled ( e.g. NaD, Ti01,2). 

The results of this analysis are shown in Figs. and |^. 
Fig. |l^ shows the run of the physical parameters, expressed as dif- 
ferentials from the reference point. In this figure the solid line indi- 
cates the full sample while triangles and diamonds represent high 
and low density samples, respectively. Fig. |lj compares the values 
of the model indices with the observed values, for the full sample. 
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Figure 9. The linearity of our adopted models for the selected indices. Index variations versus , left: metallicity (log{Z/Z..j)) at a fixed age of 10 Gyr and fixed 
[a/Fe] and [C/H], middle: log(age), and right: [a/Fe]. 
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Figure 10. Age, metallicity, a-enhancement and carbon enhancement variations as a function of velocity dispersion. A linear regression analysis has been 
carried out simultaneously on the HjS, H5, Mgl, Mg2, Mgb, Fe4383, Fe4531, Fe5270, Fe5335, G4300 and C4668 indices. The solid line represents the entire 
sample, diamonds only those objects in low density environments (l/rs < 0.5) and triangles only those in high density environments (l/r^ > 1.5). Values are 
differences with respect to those of the entire sample at a velocity dispersion of 200 km s" ' . The number of objects in each bin in CT is shown in Fig. jll] 
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Figure 12. A posteriori comparison of the models used in tlie regression analysis of Fig. |l^ witli the indices included in the analysis. The lines are the 
simultaneous solution to the whole set of indices. 



The model indices have been obtained by adding the model solu- 
tions found via equation to the reference indices. In other words 
the lines shown in Figure]! 2| represent the simultaneous solution to 
the whole set of indices. 

We note immediately that only the age is influenced by the 
density of environment; metallicity, [a/Fe] and [C/H] show no de- 
pendence on environment. We find the following trends with veloc- 
ity dispersion. 

For the sample as a whole (solid line in Fig. |^ the age rises 
significantly from -0.25 dex at 100 kms^' to a maximum near 
260 kms^'. Neither Kuntschner et al. (2001) nor Annibali et al. 
(2005) find any dependence of age on \og{a) whereas Thomas et 
al. (2005) find a significant positive slope (~0.24-0.32). Beyond 
260 kms^' there is then a slight fall to -0.07 dex in the last bin at 
320 kms~', significant at about the 2(J level. This decrease is more 
significant for galaxies in dense environments. If this decrease is 
real it may be evidence that the most massive galaxies in dense 
environments have undergone some recent star formation. 

In the range 140 — 260 kms^' galaxies in low density envi- 
ronments are ~ 0.07 dex (~ 1.5 Gyr for an elliptical of 10 Gyr) 
younger than their counterparts in dense environments for a given 



(7. Although the significance of this difference is low in any single 
bin in velocity dispersion the consistency of the difference is evi- 
dence that it is a real difference. This result is almost independent 
of the models adopted and, since the comparison is performed at 
fixed (7, it is also unlikely to be affected by the velocity dispersion 
correction or the reference to the Lick system. 

The metallicity shows a very significant, monotonic rise from 
-0.1 dex at 100 kms"' to -1-0.25 at320kms"' described weflby. 



d log(Z) 
d log(a) 



:0.76 



(6) 



and independent of environment. This is in good agreement with 
the value found by Kuntschner et al. (2001). It is 20% larger than 
that derived by Annibali et al. (2005), 30% larger than that found 
by Nelan et al. (2005) and 40% larger than that of Thomas et al. 
(2005). Although in this latter case the models show a zero-point 
offset in metallicity with respect to those used in the present study 
this should not result in a different gradient in the differential anal- 
ysis carried out here. It is more likely that the difference is due to 
differences in aperture or velocity dispersion corrections. Bernardi 
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et al. (2006) find values in the range 0.38 - 0.58, again somewhiat 
lower than our value. 

The a-enhancement rises from -0.15 dex to +0.2 dex from 100 
to 320 kms^' . We fit the linear relation; 



:0.74 



(7) 



d [a/Fe] 
d log(<7) 

This slope is about 1.5 times that found by Annibali et al. (2005) 
and about 2.5 times larger than that found by Thomas et al. (2005) 
and Kuntschner et al. (2001). It is also a factor of 2 steeper than 
that found by Bernardi et al (2006) and Nelan et al. (2005). This 
chemical enrichment pattern shows that the star formation process 
in low mass halos has occurred over a longer time scale, showing 
that the effect of feedback must be important in the baryon assem- 
bly process. Whatever the details of this feedback, environment has 
no role to play. 

The carbon enhancement instead shows a steady fall from 
+0.05 dex to -0.075 dex from the lowest to the highest velocity 
dispersion: 



-0.31 



(8) 



d [C/H] 
d log(cT) 

For more massive galaxies carbon is under-produced relative to 
other metals. This is the first time that such a trend can be quan- 
tified. 

We see from Fig. |^ that despite the trends seen with velocity 
dispersion in Fig, there is nonetheless a large spread in the index 
values for any given value of o. Therefore there exist galaxies with 
CJ~ lOOkms"' that are as old and metal rich as as those with CT ~ 
300 kms~' . However, the fact that the range of Hj8 and H5 values 
is much smaller at large velocity dispersion than at low velocity 
dispersion is evidence that the converse is not true: galaxies with 
(T ~ 300 kms^' are never as young and metal poor as the typical 
galaxy with a ~ 100 km s^'. 



5 DISCUSSION 

The current sample of early-type galaxies has been selected in such 
a way as to avoid the redshift/density bias that is introduced by tra- 
ditional selection techniques when applied to a survey with a lim- 
ited sky coverage. With a sample so defined we find that there is no 
difference in either the redshift distribution or the z-band luminos- 
ity distribution of early-type galaxies in more or less dense envi- 
ronments. We do however find that environment affects the distri- 
bution of velocity dispersion, with the mean <J rising with the local 
density of the environment. This permits the analysis of the sam- 
ple as a whole without considering separate bins of either redshift 
or galaxy luminosity, and the consideration only of the distribu- 
tion of CT within restricted ranges. Furthermore the large number 
of galaxies allows a robust estimate of the mean values of the in- 
dices at increasing velocity dispersion, even when sub-samples of 
objects belonging to different environments are selected. In fact the 
presence of environmental effects has already been qualitatively in- 
ferred from the statistical analysis described in section 12 



From the quantitative analysis performed with new models, in 
which we have included all the major dependencies from the phys- 
ical parameters, we infer that both the enhancement of a-elements 
and the global metallicity increase with galactic velocity disper- 
sion. Instead, contrary to previous studies, we find that the age 
initially increases but then flattens, and possibly then falls at the 
highest velocity dispersions. 



While some of these trends have already been uncovered by 
previous studies (Bressan et al. 1996, Thomas et al., 2005) it is 
only by means of a careful statistical analysis of a large data sam- 
ple and its interpretation with the most up-to-date SSPs presented 
here, that it is possible to obtain perhaps the cleanest scenario so far 
for the formation and evolution of early type galaxies in different 
environments. Despite broad agreement with recent studies based 
on similarly large samples (Nelan et al., 2005, Bernardi et al., 2006, 
Smith et al., 2006) one key difference emerges; that is our finding 
that the only environment sensitive parameter is age. The above 
authors also find slightly enhanced a abundances in denser envi- 
ronments. Whether this difference is a result of sample selection 
criteria or diverse approaches to data analysis is not clear. 

The general picture which emerges from our work is that at 
fixed mass (cr), galaxies formed first in dense environments (that 
now correspond to clusters) and ~ 1.5 Gyr later in the less dense 
field, but with a similar star formation time-scale. If the early col- 
lapse in high density peaks happened at a redshift of say z ^ 6 then, 
a 1.5 Gyr delay is compatible with a late collapse at redshift of 
z — 2.5. Lower mass galaxies (a ~ 100 kms^') formed ~ 4 Gyr 
more recently in both environments, which would correspond to 
a mean formation redshift of ~ 1.1 if those with o ~ 220 kms^' 
formed at z = 6. The star formation rate that characterised the mass 
assembly was more rapid for more massive objects (from the in- 
crease in [oj/Fe] with cr). 

The observations are at odds with the predictions of usual hier- 
archical models of galaxy formation, where smaller halos typically 
collapse earlier than the larger ones (Lacey & Cole, 1994), indicat- 
ing that an important ingredient is still missing in these models. If 
we assume a standard initial mass function (IMF) i.e. not biased 
toward high mass stars, the trend of a-enhancement with velocity 
dispersion suggests that the stars that make up more massive early- 
type galaxies formed over characteristically shorter timescales than 
those that make up lower mass objects. Furthermore, in the absence 
of an effect of environment on the a-enhancement, we are led to 
conclude that the timescale of a typical star formation episode is 
not affected directly by the environment, only by the mass of the 
object in question. Since at constant efficiency the metal enrich- 
ment is larger when the star formation process lasts longer, this 
implies that the baryon assembly must have been much less effi- 
cient in smaller systems. This is strong evidence for the importance 
of the energetic feedback on the star formation process, as recently 
shown by Granato et al. (2004). In these models the star formation 
process in smaller systems is rendered less efficient by the super- 
nova feedback with the result that, though starting earlier, the star 
formation lasts longer, giving rise to a lower metal enrichment. 

Although slower star formation with stronger feedback in low- 
mass systems explains the observed trends in metallicity and a- 
enhancement it does not get over the problem of how most of the 
star formation occurred after the mass was assembled. If this were 
not the case we would not observe a rising [a/Fe] with velocity 
dispersion. Given that gas-rich mergers are expected to be accom- 
panied by bursts of star formation it is difficult to see how the star 
formation can be delayed if multiple-mergers are responsible for 
the mass assembly. Conversely, a small number of mergers would 
produce systems with a significant net angular momentum. 

It is not possible from the present analysis to quantify the ef- 
fects of major mergers, though their efficiency should be higher 
in the kind of low-velocity encounter typical of the field than in 
the high-velocity encounters expected in clusters. It could then also 
be possible that the formation process in field early type galaxies 
is modulated by the occurrence of major mergers. Annibali et al. 
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(2005) have shown that in that case the fraction of stars formed 
during the burst is less than 25%. 

Finally, we find that [C/H] decreases at increasing metal con- 
tent. This carbon anomaly is suggested by comparison of observed 
index values (Fig. ^ and is thus independent of the models. 

Models in which [C/H] scales with the metal content do not 
provide consistent results between indices with different sensitiv- 
ity to this element. The variation of [C/H] between low and high a 
galaxies is approximately 0.125 dex, independent of environment. 
It is worth noting that several indices, usually adopted for this kind 
of theoretical analysis, have a more or less pronounced dependence 
on the carbon abundance (including Mgl, Mg2 and to a lesser ex- 
tent Mgb, see e.g. Korn et al. 2005). It is thus important to perform 
the analysis with suitable models that account for this effect. This 
abundance trend may suggest a lower production of carbon in more 
metal rich, high mass stars. But it may also indicate a lower C yield 
from intermediate mass stars as the 'hot bottom burning' process 
becomes more efficient at increasing metal content (Marigo, Bres- 
san & Chiosi 1998) 



6 CONCLUSIONS 

We have carried out a detailed analysis of the spectral line in- 
dices of a sample of 3614 early-type galaxies selected form the 
Sloan Digital Sky Survey. Our selection technique avoids the red- 
shift/environment bias introduced by traditional selection criteria 
when applied to a survey with partial sky coverage. We have quan- 
tified the density of environment for each object and have in- 
vestigated the effects of environment on age, metallicity and a- 
enhancement via a regression analysis of the whole sample. Be- 
cause the carbon abundance has a strong effect on the strength of 
many indices we have used revised models which take account also 
of the carbon abundance explicitly. Thus, our models take into ac- 
count the main parameters that drive the set of indices adopted in 
the analysis; age, metallicity, a-enhancement and carbon content. 

Our results for early-type galaxies can be summarised as fol- 
lows. 

• Neither the z-band luminosity nor the redshift distribution are 
a function of galaxy environment out to a redshift of 0.1. 

• The mean velocity dispersion is higher in denser environ- 
ments. 

• The only effect of enviroimient is on the mean age of early- 
type galaxies for a given velocity dispersion. Environment has no 
effect on metallicity, a-enhancement or [C/H]. 

• The star formation episodes responsible for the assembly of 
the baryonic mass of early-type galaxies were of shorter duration 
in objects which have a higher present day mass. Environment did 
not influence the duration of these episodes. 

• The youngest objects are those with the lowest masses, both 
in high and low density environments. 

• Objects in the field with all but the highest velocity dispersion 
are, on average, 1.5 Gyr younger than their cluster counterparts. 
The same trend is also obtained by adopting the models by Thomas 
et al. (2005). 

• Carbon becomes increasingly under-abundant in higher mass 
objects. 

The trend we find of increasing a-enhancement with galaxy 
mass is consistent with a pure monolithic collapse model. However, 
the more prolonged star formation in lower mass objects should 
also produce a higher metallicity in these objects, which is not seen. 



A pure monolithic collapse model is therefore rejected. In a simple 
hierarchical merger scenario the a-enhancement trend is not repro- 
duced. The a-enhancement trend could be explained in a hierar- 
chical model if small mass halos which formed within a large dark 
matter halo formed their stars more rapidly than equivalent mass 
objects in a smaller parent halo. However, such a picture would 
predict a dependency of the star formation timescale (and thus a- 
enhancement) also on environment. We see no such effect with en- 
vironment. Galaxy formation in which the star formation process 
is dependent only on the mass of the galaxy is not consistent with 
either simple monolithic collapse or hierarchical merging. 

A rising a-enhancement and metallicity with velocity disper- 
sion requires a feedback mechanism that results in, 1) a slower star 
formation in low mass objects, and 2) a less efficient star formation 
in low mass objects. The latter ensures that the more prolonged 
star formation in low mass objects does not result in more metals. 
A model with mass-dependent feedback has been investigated by 
Granato et al. (2004) 

Small mass halos then, can collapse before larger mass ha- 
los, only if the star formation within them is inhibited by such a 
feedback mechanism, such that, in the absence of mergers, their 
star formation proceeds at an extremely slow rate. When mergers 
take place the increased mass reduces the feedback and the rate 
and efficiency of star formation increases. The greater the merger 
mass the more the star formation rate is increased and the sooner it 
ends. The result is that the mean age of the stars formed increases 
with mass even though the merger epoch can be earlier for lower 
mass systems. There is therfore an important distinction between 
the epoch of assembly of the dark halo and the mean age of the 
stars in a galaxy which we measure spectroscopically. The differ- 
ence is largest for low mass galaxies. In environments with a high 
density the merger process occurred --^1.5 Gyr earlier than in low 
density environments and this timing seems to be the only influence 
of the environment. 
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APPENDIX A: ROBUSTNESS OF RESULTS TO 
VELOCITY DISPERSION CORRECTION OF INDICES 

Our findings, that the index values for early-type galaxies are af- 
fected by environment, are not due to an error in our velocity dis- 
persion correction. If an error in velocity dispersion correction in- 
troduced an erroneous gradient component into the index- CT rela- 
tions this would not alter the separation between the 2 lines in Rig. ^ 
because the Monte-Carlo simulations use the measured index-c re- 
lations in order to establish how values of a are translated into val- 
ues of the index. In order to show this we carried out two tests. In 
the first we introduced an error term into the gradient of the index- 
es relations and re-computed the simulations illustrated in fig. |^ 
No differences were seen that extended beyond the error bars in 
this figure. Secondly, we constrained the CT-distributions of the sub- 
samples in each environment to be identical to that of the compar- 
ison sub-sample ( 1 /r^ < 0.5) by randomly selecting objects until 
the number of objects per bin was in constant proportion to the 
numbers in the comparison sample. This of course leads to the ex- 
clusion of many objects but as this is at random only the size of the 
errors on individual points will be affected. Reproducing Rig. |8| in 
this way also did not change the index difference between data and 
simulations. The results of this test are shown in Rig. where we 
have subtracted the simulation from the data to illustrate the effect 
of environment alone on the line indices. In short, our results on 
the effect of environment are robust to errors in velocity dispersion 
correction because it is not dependent on the form of the index-cT 
relation. 

In order to test whether these differences affect only the low 
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Figure Al. Effect of environment on relative line index strength free of the effect of the index-CT relation. Similar to Fig |^ but shown is the difference between 
the data and simulation. The solid line is for the the sub-samples exactly as defined in Fig ^ and the dotted line is for sub-samples which have been constrained 
(by random sampling in bins of a) to have similar velocity dispersion distributions. The error bars are not shown for the second line but are similar to those 
of the solid line. The similarity of the lines shows that neither an error in velocity dispersion correction nor the differences in velocity dispersion distribution 
between sub- samples can explain the effect we find with environment. 
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mass objects or are general to the early-type population as a whole 
(at least for objects brighter than -20.45 in R-band) we have re- 
peated our analysis for 3 bins in velocity dispersion chosen to rep- 
resent low, intermediate and high mass systems: CT < 160 kms^', 
160 ^ (7 ^ 230 kms^' and <J > 230 kms^'. The luminosity and 
redshift distributions of these sub-samples are shown in Fig. A2, 



where it can be seen that they are similarly distributed. The results 
of the analysis for these 3 ranges of velocity dispersion are shown 
in the rightmost 3 panels for each index of Fig. H. 

Bearing in mind the larger errors on the results for the sam- 
ple when divided by velocity dispersion (especially for the low 
and high sigma sub-samples) we find that most of the trends seen 
for the sample as a whole are also seen for the sub-sample with 
a < 160 kms^' . Indeed, the magnitude of the environmental index 
differences is larger for these low-mass objects for all the indices 
in which some effect is seen for the sample as a whole. The sub- 
sample with tr > 230 kms^' shows effects with environment that 
are generally weaker than those for the low-mass sub-sample and 
for the sample as a whole such that the magnitude of the environ- 
mental effect falls from panel 2 to panel 4 of Fig. ^. CNl and Mgl 
show this trend especially clearly. However there are exceptions to 
this trend. G4300, Hb, Fe5015 and Fe5270 seem to show stronger 
trends with environment for more massive objects although none of 
the trends are as convincing as those in the opposite sense. 
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Figure A2. Distribution of redshift for vaiious density environments for 3 ranges of the velocity dispersion. The lines and labelling are as in fig. ^ Left: 
a < 160 kms"', middle: 160 ^ ff ^ 230 kms"', right: a > 230 kms"'. The curves are normalised to the number of objects in the low density comparison 
sample which contains, 185, 473, and 68 objects respectively for the 3 ranges. 
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Figure A3. As for fig. for the distribution of luminosity. 



